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Abstract 

We discuss reflection factors for purely elastic scattering theories and relate 
them to perturbations of specific conformal boundary conditions, using re¬ 
cent results on exact off-critical ^-functions. For the non-unitary cases, we 
support our conjectures using a relationship with quantum group reductions 
of the sine-Gordon model. Our results imply the existence of a variety of 
new flows between conformal boundary conditions, some of them driven by 
boundary-changing operators. 


1 Introduction 


There are two general approaches to the study of integrable boundary quantum 
field theories: infrared and ultraviolet. From the infrared point of view, solving 
the boundary Yang-Baxter and bootstrap constraints allows sets of reflection 
factors to be associated with various bulk scattering theories. However, there are 
many ambiguities left by these constraints, which cannot be lifted without the 
introduction of more selective criteria. This is to be expected, since a given bulk 
theory may admit many different integrable boundary conditions. 

Alternatively, a theory might be defined via a classical Lagrangian, or as a 
perturbation of a boundary conformal field theory. In either case, it is natural 
to ask for a ‘UV/IR dictionary’, giving the relationship between the infrared and 
ultraviolet specifications of the model. In particular, this would reveal which 
among the infinite families of solutions to the boundary Yang-Baxter and boot¬ 
strap equations are physically realised, and to which boundary conditions they 
correspond. 

Progress on this issue includes the relationship between the parameters in 
the boundary reffection factors and the couplings in the boundary sine-Gordon 
Lagrangian derived by Aliosha Zamolodchikov [1] (checked using the boundary 
TCSA [2] in [3]), and various perturbative results for cases where the theory ad¬ 
mits a Lagrangian definition (see for example [4, 5]). However, these methods 
cannot be readily adapted to the study of general perturbed boundary confor¬ 
mal field theories. In particular, a boundary analogue of the thermodynamic 
Bethe ansatz (TBA) c-function, which allows bulk S-matrices to be identified 
with specific perturbed conformal field theories [ 6 , 7], has been lacking. An ob¬ 
vious candidate is the 51 -function defined by Affleck and Ludwig in [ 8 ], but the 
initial proposal for a TBA-like equation for a fully off-critical version of g, taking 
as input just the bulk S-matrix and the boundary reflection factors [9], turned 
out to be incorrect [ 10 ]. 

Recently it has been shown that the situation can be remedied [11]. The de¬ 
tailed work in [11] concentrated on the boundary scaling Lee-Yang model, for 
which the relationship between microscopic boundary conditions and reflection 
factors had already been established [2]. In this paper we extend the investiga¬ 
tions of [11] to a collection of theories for which boundary UV/IR relations have 
yet to be found, namely the minimal purely elastic scattering theories associated 
with the ADET series of diagrams [12, 7, 13, 14, 15, 16, 17]. The bulk S-matrices 
of these models have long been known, but less progress has been made in as¬ 
sociating solutions of the boundary bootstrap equations with specific perturbed 
boundary conditions. We present a collection of minimal reffection factors for 
the ADET theories, and test them by checking the ^f-function ffows that they 
imply. We also show how these reflection factors can be modified to incorporate 
a free parameter, generalising a structure previously observed in the Lee-Yang 
model [2]. This enables us to predict a number of new ffows between conformal 
boundary conditions. 
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In more detail, the rest of the paper is organised as follows. Bulk and boundary 
ADET scattering theories are discussed in sections 2 and 3, with section 3 con¬ 
taining a full set of reflection factors fulfilling certain minimality criteria. These 
solutions have no free parameters, and in section 4 the minimality requirement 
is relaxed and non-minimal one-parameter families of solutions are proposed to 
describe simultaneous bulk and boundary perturbations of boundary conformal 
field theories by relevant operators. In section 5 the one-parameter families of 
reflection amplitudes for the models are alternatively deduced by considering 
special reductions of the boundary sine-Gordon model. 

After a discussion of the exact off-critical < 7 -function in section 6 , the phys¬ 
ical consistency of the minimal solutions is checked in section 7, by comparing 
the predictions for the ultraviolet values of the g'-function obtained from the re¬ 
flection factors with conformal field theory results. An additional check in first 
order perturbation theory is performed for the 3-state Potts model in section 8 . 
Section 9 summarises the ^r-function predictions for the one-parameter families of 
reflection factors and discusses their UV/IR relations. Our conclusions are given 
in section 10 , and an appendix gathers together various results concerning the 
fields and function values in the diagonal 0 i x 01/02 coset models which are 
needed in the main text. 

2 The ADET family of purely elastic scattering theo¬ 
ries 

Purely elastic scattering theories are characterised by the property that their S- 
matrices are diagonal. The scattering of particles a and h with relative rapidity 
9 is then described by a single function Sab{9), which is a pure phase for physical 
rapidities. The unitarity and crossing symmetry conditions simplify to 


Sab{0)Sab{-e) = 1, 

(2.1) 

Sab{9) = 5,-,(i7r-0) 

(2.2) 


respectively, where b is the antiparticle of b. The Yang-Baxter equation is trivially 
satished, but, as stressed by Zamolodchikov [18, 19], the bootstrap still provides a 
useful constraint. This states that whenever an on-shell three-point coupling 
is nonzero, the S-matrix elements of particles a, b and c with any other particle 
d satisfy 

SdM = SdaiO - iUL)Sdb{0 + (2.3) 

where the ‘fusing angles’ are related to the masses of particles a, b and c by 

m2 = m 2 + Mfe2 + 2MaMb cos (2.4) 

and U = TT — U. These angles also appear as the locations of certain odd-order 
poles in the S-matrix elements. More details of the workings of the bootstrap can 
be found in [18, 19] and, for example, [20]. 
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The purely elastic scattering theories that we shall treat in this paper fall 
into two classes. The first class associates an S-matrix to each simply-laced Lie 
algebra 0 , of type A, D or E [12, 7, 13, 14], These S-matrices are minimal, in 
that they have no zeros on the physical strip 0 < < tt, and one-particle 

unitary, in that all on-shell three-point couplings, as inferred from the residues of 
forward-channel S-matrix poles, are real. They describe particle scattering in the 
perturbations of the coset conformal field theories 0 i x 0 i /02 by their ( 1 , l,ad) 
operators, where g is the affine algebra associated with g. The unperturbed 
theories have central charge 


c = 


(h + 2) ’ 


(2.5) 


where rg is the rank of g, and h is its Coxeter number. These UV central charges 
can be recovered directly from the S-matrices, using the thermodynamic Bethe 
ansatz [ 6 , 7]. 

The ADE S-matrices describe the diagonal scattering of rg particle types, 
whose masses together form the Perron-Frobenius eigenvector of the Cartan ma¬ 
trix of 0 . This allows the particles to be attached to the nodes of the Dynkin 
diagram of g. Each S-matrix element can be conveniently written as a product 
of elementary blocks [ 12 ] 


sink (- -I- 

W = (X-1)U+1). W(.) = -A|^ 


as 


Sab = n ’ 


( 2 . 6 ) 

(2.7) 


for some index set Aab- Note that 


( 0 ) = 1 , {h) = -l, (-x) = (x) \ (x± 2 /i) = (x). 


( 2 . 8 ) 


The notation (2.6) has been so arranged that the numbers x are all integers. The 
sets Aab are tabulated in [ 12 ]; a universal formula expressing them in geometrical 
terms was found in [14], and is further discussed in [15, 16]. 

The S-matrices of the second class [ 21 , 22 , 23] are labelled by extending the set 
of ADE Dynkin diagrams to include the ‘tadpole’ T^. They encode the diagonal 
scattering of r particle types, and can be written in terms of the blocks ( 2 . 6 ) with 
h = 2r+l [24]; 

Sab= n (2.9) 

l=\a—b\-\-l 
step 2 

These S-matrix elements are again minimal, but they are not one-particle unitary, 
reflecting the fact that they describe perturbations of the non-unitary minimal 
models M 2 , 2 r-\- 3 , with central charge c = —2r(6r -|- 5)/(2r -|- 3). The perturbing 
operator this time is (pis . The Tr S-matrices are quantum group reductions of the 


3 



sine-Gordon model at coupling 0^ = 167r/(2r+3) [21], a fact that will be relevant 
later. 

A self-contained classification of minimal purely elastic S-matrices is still lack¬ 
ing, but the results of [17] single out the ADET theories as the only examples 
having TEA systems for which all pseudoenergies remain finite in the ultraviolet 
limit. 

The ADET diagrams are shown in figure 1, with nodes giving our conventions 
for labelling the particles in each theory. Eor the Dj. theories, particles r—1 and 
r are sometimes labelled s and s', or s and s, for r even or odd respectively. 


• — • — • .... • — • 

1 2 3 r-1 r 

Fig. la: Ar 



Fig. lb: Dr 


• • • • ■ • • #0 

12 3 r 

Fig. Ic: Tr 
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Fig. Id: Eq 
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Fig. le: E'j 



3 Minimal reflection factors for purely elastic scatter¬ 
ing theories 

The scattering of particles by a boundary in an integrable quantum field theory is 
described by a set of reflection factors. In this paper we shall restrict ourselves to 
situations where not only the bulk, but also the boundary scattering amplitudes 
are purely elastic. The reflection factors are then rapidity-dependent functions 
R0O), one for each particle type a in the theory. Unitarity [25, 26] and crossing- 
unitarity [26] constrain these functions as 

R09)Rai-O) = 1, (3.1) 

R09)R09-i7r) = Saa{29). (3.2) 

In addition, whenever a bulk three-point coupling is nonzero, a boundary 
bootstrap equation holds [25, 26]: 

R09) = R09 + iUl)R09 - iU^0Sa029 + iUl - iU0) . (3.3) 

For a given bulk S-matrix there are infinitely-many distinct sets of reflection 
factors consistent with these constraints, because any solution can be multiplied 
by a solution of the bulk bootstrap, unitarity and crossing equations to yield 
another solution [27]. To identify a set of reflection factors as physically relevant 
some more information is needed. A common basis for the conjecturing of bulk 
scattering amplitudes is a ‘minimality hypothesis’, that in the absence of other 
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requirements one should look for solutions of the constraints with the smallest 
possible number of poles and zeros. In this section we show how the same principle 
can be used to find sets of boundary amplitudes, one for each purely elastic S- 
matrix of type A, D or E. The amplitudes given below were at first conjectured [28] 
as a natural generalisation of the minimal versions of the Ar affine Toda field 
theory amplitudes found in [29, 30]. The reasoning behind these conjectures will 
be explained shortly. More recently, Fateev [31] proposed a set of reflection factors 
for the affine Toda field theories. These were in integral form, and not all matrix 
elements were given. However, modulo some overall signs and small typos, the 
coupling-independent parts of Fateev’s conjectures match ours. (This has also 
been found by Zambon [32], who obtained equivalent formulae to those recorded 
below taking Fateev’s integral formulae as a starting point.) 

Unitarity and crossing-unitarity together imply that the reflection factors 
must be 27ri periodic; unitarity then requires that they be products of the blocks 
(x) introduced in the last section. The crossing-unitarity constraint is then key 
for the analysis of minimality. Each pole or zero of Saa{‘^Q) on the right hand 
side of (3.2) must be present in one or other of the factors on the left hand side 
of that equation. This sets a lower bound on the number of poles and zeros for 
the reflection factor Ra{0). 

To exploit this observation, it is convenient to work at the level of the larger 
blocks {x}, the basic units for the bulk bootstrap equations [12, 14]. We define 
two complementary ‘square roots’ of these blocks as 

(a:) = (£^) , (J) = (s+1) (3,4) 

and record their basic properties 

{x){9) = {x){9 + ill) (3.5) 

and 

{x){9){x){9) = {x]{29). (3.6) 

The crossing-unitarity equation can then be solved, in a minimal fashion, by any 
product 

Ra= W L (3.7) 

where each factor fx can be freely chosen to be (x) or (x), modulo one subtlety: 
since (0) = 1, minimality requires that any factor /i be taken to be (1) rather 
than (1). In fact there is always exactly one such factor for the diagonal S-matrix 
elements relevant here [15].* 

*We could also set {x) = (^) {2^), {x) = (24^^)“^ throughout, but this 

would break the pattern previously seen for the Ar theories, and turns out not to fit with the 
(/-function calculations later. 
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There remain the bootstrap equations (3.3). To treat these we need to know 
how the blocks {x){6) and {x){6) transform under general shifts in 9. As in 
[14, 15], these shifts are best discussed by defining 

(x)+ (9) = sinh (^- + (3.8) 

and then 

W+ = W)+(^^)i" . (J>+ = (T)+(2^)i" (3-9) 

so that (x) = (x)+/(—x)+, (x) = (x)+/(—x)+ and (x) = (x)+/(—x)+. Once a 
putative set of reflection factors (3.7) has been decomposed into these blocks, it is 
straightforward to implement the boundary bootstrap equations (3.3) using the 
properties 

{x)+{9 + my/h) = (x+y) + (6') , {x)+{9 + my/h) = {x+y)+{9) . (3.10) 

For the ADE theories it turns out that the boundary bootstrap equations can 
all be satisfied by minimal conjectures of the form (3.7), and that the choice of 
blocks is then fixed uniquely. (From the point of view of the bootstrap equations 
alone, an overall swap between (x) and (x) is possible, but the requirement that 
/i = (1) fixes even this ambiguity.) The final answers are recorded below. 

The reflection factors for this case were given in [29]. In the current notation they 
are 

2a-1 

i?, = n (^) (3-11) 

X = 1 

X odd 


Dr 

The reflection factors, for a = 1,..., r — 2 are 

2a—1 2a—1 2a—5 

^a = n n {h — x) (/i — X — 2 ), for a odd 

x=l X = 1 X = 1 

X odd step 4 step 4 

2a—1 2a—3 ^ ^ 

(x) {h — x){h — X — 2), for a even 


X=1 X=1 

X odd step 4 


(3.12) 
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while for a = r — 1 and r we have 

2r-5 

Ry-i = Ry. = {x ), for r odd (3.13) 

step 4 
2r-3 

Rr-i = (^)) even. 

31 = 1 
step 4 

E, 

Ri 
R 2 
Rs 

Ri 
E-r 

Ri 
R 2 
R 3 

Ri 

R 5 
Rq 


Es 

Ri = (1)(11)(19)(29)^ 

R^ = (1)(7)(11)(13)(27)(19)(23)(29)^ 

Rii = (1)(3)(9)(11)(11)(^13)(17)(19)(29)(21)(27)(^29) 

Ri = (l)(5)(7)(9)(llKll)(13)a5)(15)a7)(19)(19)(21)(23)(25)(29) 

R^ = (1)(3)^5)(7)(9K9)(11)2(11)(13)(13)(15)2(17)(17)(19)2(19)(21)(21) 

(23)(25)(27)(29) (3.16) 

Ra = (I)(3)(5)(7)j7)(9)2(n)2(ri)(13)2(r3)(15)(r5)(17)2(r7)(19)(r9)2 
(21)2(23)(23)(25)(27)(29) 

= (I)(3)(5)2(7)2(7)(9)2(9)(ll)2(ri)2(13)^(f3)(15)2(r5)2(17)^(r7)(19)2 
(19)2 (21)2 (2~1) (23) (2'3)2 (25)2 (2~7) (29) 

= (I)(3)2(5)2(5)(7)^(7)(9)^(9)2(ll)^(ri)2(l3)3(r3)3(15)^(r5)2(17)^(r7)^ 
(19)^(f9)2(21)2(2l)3(23)^(2'3)(25)(2'5)2(27)2(f9) 


= (1)(9)(17) 

= (1)(7)(11)(17^ 

= (1)(5)(7)(9)(11)(13)(17)^ 

= (1)(3)(7)(9)(9)(11)(15)(1J) ^ (3.15) 

= (1)(3)(5)(7)(7K9)2(11)(11)(1J1)(15)(17) 

= (l)(3)(5)2(7K7)(9)2(9)(llKll)(13)2jl5)(17^ 

= (1)(3)2(5)2(5)(7)3(7)(9)2(9)2(11)3(11)(13)(13)2(15)2(17) 


= R-i = (1)(7) 

= (1)(5)(7)(11) ^ (3.14) 

= i?3 = (1)(3K5)(7)^7)(9)^ 

= (1)(3)2(5)2(5)(7)2(7)(9)(9)(11) 
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For the T series the story is different: it is not possible to satisfy the bound¬ 
ary bootstrap equations with a conjecture of the form (3.7). This means that the 
minimal reflection factors for these models are forced by the bootstrap to have 
extra poles and zeros beyond those required by crossing-unitarity alone. Our gen¬ 
eral proposal will be given in eq. (5.4) below, but the situation can be understood 
using the boundary Ti, or Lee-Yang, model: the minimal reflection factor found 
in [2] for the single particle in this model bouncing off the |1) boundary is 

= (5)(i)(|)“' • (3-17) 

The simpler function (^)(|)~^ would have been enough to satisfy crossing-unitarity, 
but then the boundary bootstrap would not have held, and so (3.17) really is a 
minimal solution. This observation fits nicely with the ^(-function calculations to 
be reported later: had the minimal reflection factors for the theories fallen into 
the pattern seen for other models, there would have been a mismatch between 
the predicted UV values of the gf-functions and the known values from conformal 
field theory. 

4 One-parameter families of reflection factors 

The minimal reflection factors introduced in the last section have no free param¬ 
eters. However, combined perturbations of a boundary conformal field theory by 
relevant bulk and boundary operators involve a dimensionless quantity - the ratio 
of the induced bulk and boundary scales - on which the reflection factors would 
be expected to depend. To describe such situations, we must drop the minimality 
hypothesis, and extend our conjectures. 

A first observation, rephrasing that of [27], is that given any two solutions 
Ra{0) and R'a{d) of the reflection unitarity, crossing-unitarity and bootstrap rela¬ 
tions (3.1), (3.2) and (3.3), their ratios Za{9) = Ra{0) / R^iS) automatically solve 
one-index versions of the bulk unitarity, crossing and bootstrap equations (2.1), 
(2.2) and (2.3): 


Za{e)Za{-e) = 1 , Za{e) = Z-aiiTT - 9) , (4.1) 

Z,{9) = Za{9 - iUl)Z,{9 + iU^,) . (4.2) 

The minimal reflection factors Ra{9) can therefore be used as multiplicative ‘seeds’ 
for more general conjectures R'a{9) = Ra{9) / Za{9), with the parameter-dependent 
parts Za{9)~^ constrained via (4.1) and (4.2). An immediate solution is z'^f^{9) = 
Sda{9) for any (fixed) particle type d in the theory, where we have used the ket 
symbol to indicate that the label \d) might ultimately refer to one of the possible 
boundary states of the theory. However, this does not yet introduce a parameter. 
Noting that a symmetrical shift in 9 preserves all the relevant equations, one 
possibility is to take 


z[^^^\9) = Sda{9 + C)Sda{9 - C), 


(4.3) 



with C at this stage arbitrary. This is indeed the solution adopted by the bound¬ 
ary scaling Lee-Yang example studied in [2], This model is the r = 1 member of 
the Tj. series described earlier, and corresponds to the perturbation of the non¬ 
unitary minimal model A ^25 by its only relevant bulk operator, ip, of conformal 
dimensions The minimal model has two conformally-invariant 

boundary conditions which were labelled |1) and |$) in [2], The |1) boundary 
has no relevant boundary fields, and has a minimal reflection factor, given by 
(3.17) above. On the other hand, the |<h) boundary has one relevant boundary 
field, denoted by 4 >, and gives rise to a one-parameter family of reflection factors 

R\^\e) = E}^\e)/z\^\9) (4.4) 

where the factor (9) has exactly the form mentioned above: 

Z\’’\9) = S{9 + if (6+3))5(0 - if (6+3)) (4.5) 

where S{9) is the bulk S-matrix, and the parameter 6 can be related to the 
dimensionless ratio +^/A of the bulk and boundary couplings A and fi [ 2 , 10]. 
(Since there is only one particle type in the Lee-Yang model, the indices a, d and 
so on are omitted. We also changed the notation slightly from that of [2] to avoid 
confusing the parameter 6 with a particle label.) ^ 

As an initial attempt to extend (4.4) to the remaining ADET theories one 
could therefore try 

= (4.6) 

with Z^a’^\0) as in (4.3). This manoeuvre certainly generates mathematically- 
consistent sets of reflection amplitudes, but in the more general cases it is not the 
most economical choice. Consider instead the functions obtained by replacing the 
blocks {x} in (2.7) by the simpler blocks (x) [33]: 

= n (^) • (4-7) 

For the Lee-Yang model, S^{9) coincides with S{9), but for other theories it has 
fewer poles and zeros. The bootstrap constraints are only satisfied up to signs, 
but if is used to define a function as 

= sl{9 + - i^c) (4.8) 

then these signs cancel, and thus (4.8) provides a more “minimal” generalization 
of the family of Lee-Yang reflection factors which nevertheless preserves all of 
its desirable properties. Setting d equal to the lightest particle in the theory 
generally gives the family with the smallest number of additional poles and zeros, 
but we shall see evidence later that all cases have a role to play. 

^There are also boundary-changing operators, but these were not considered in [2, 10]. 
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The normalisation of the shift was changed in passing from (4.3) to (4.8); this 
is convenient because, as a consequence of the property 

(x - C)(0) x(x + C)(0) = {x){e + ifC) X {x){e - ifC ), (4.9) 

we have 

zr>= n {x -C){x + C) . (4.10) 

The factors therefore coincide with the coupling-constant dependent 

parts of the affine Toda field theory S-matrices of [12, 7], with C related to the 
parameter B of [12] by C = 1 — i?. 


5 Boundary Tr theories as reductions of boundary sine- 
Gordon 


The reflection factors presented so far are only conjectures, and no evidence has 
been given linking them to any physically-realised boundary conditions. The 
best signal in this respect will come from the exact gf-function calculations to be 
reported in later sections. However, for the theories, alternative support for 
our general scheme comes from an interesting relation with the reflection factors 
of the sine-Gordon model. 

In the bulk, the theories can be found as particularly-simple quantum group 
reductions of the sine-Gordon model at certain values of the coupling, in which the 
soliton and antisoliton states are deleted leaving only the breathers [21]. Quantum 
group reduction in the presence of boundaries has yet to be fully understood, 
but the simplifications of the cases allow extra progress to be made. This 
generalises the analysis of [2] for Ti, the Lee-Yang model, but has some new 
features. 

The boundary S-matrix for the sine-Gordon solitons was found by Ghoshal 
and Zamolodchikov in [26], and extended to the breathers by Ghoshal in [34]. To 
match the notation used above for the ADE theories, we trade the sine-Gordon 
bulk coupling constant (3 for 



(5.1) 


Ghoshal and Zamolodchikov expressed their matrix solution to the boundary 
Yang-Baxter equation for the sine-Gordon model^ in terms of two parameters ^ 
and k. However for the scalar part (which is the whole reflection factor for the 
breathers) they found it more convenient to use rj and 'd, related to ^ and k by 


cos ri cosh d = —— cos ^ , cos^ rj + cosh^ d = 1 + 
k 




Glso found by de Vega and Gonzalez Ruiz [35] 


(5.2) 


10 



Ghoshal-Zamolodchikov’s reflection factor for the breather on the |r/, i?) bound¬ 
ary can then be written as 



R\y\e) = B!0\9)d^\9) 

(5.3) 

where 

p(a) (1) {a + h) {1) {1 + h) 

° (« + f) + 

is the boundary-parameter-independent part and 

(5.4) 


R^^\9) = S^^\rj,9)S^^\i^,9) 

(5.5) 

with 

a-l (2u h 

n , L, 

l=l-a V TT 2 ‘ ; 

(5.6) 


step 2 


contains the dependence on rj and i?. In these formulae we used the same blocks 
(x) as defined in (2.6) for the minimal ADE theories, with h now given by (5.1). 

Now for the quantum group reduction. In the bulk, suppose that 0^ is such 
that 

h = 2r + l, reN. (5.7) 

At these values of the coupling, Smirnov has shown [21] that a consistent scatter¬ 
ing theory can be obtained by removing all the solitonic states, leaving r breathers 
with masses Ma = • The scattering of these breathers is then given 

by the S-matrix (2.9). In order to explain all poles without recourse to the 
solitons, extra three-point couplings must be introduced, some of which are neces¬ 
sarily imaginary, consistent with the models being perturbations of nonunitary 
conformal held theories. In the presence of a boundary, these extra couplings give 
rise to extra boundary bootstrap equations, which further constrain the boundary 
reflection factors and impose a relation between the parameters t] and i). This 
can be seen by a direct analysis of the boundary bootstrap equations but it is 
more interesting to take another route, as follows. 

We first recall the observation of [2], that the reflection factors for the bound¬ 
ary Lee-Yang model match solutions of boundary Yang-Baxter equation for the 
sine-Gordon model at 

^ ^ ioo . (5-8) 

The Lee-Yang case corresponds to r = 1, /i = 3, but we shall suppose that the 
same constraint should hold more generally (see also [9]). The condition must be 
translated into the (77,7?) parametrisation. One degree of freedom can be retained 
by allowing k to tend to infinity as the limit (5.8) is taken. The relations (5.2) 
then become 

cos rj cosh 7? = A , cos^ rj + cosh^ 7? = 1. (5-9) 
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The constant A can be tuned to any value by taking ^ and k to infinity suitably, 
and so the first equation in (5.9) is not a constraint; solving the second for i-d, 


id 


TT 


r/ + - + (r-d)7r , 


d G Z. 


(5.10) 


(The shift by r is included for later convenience.) This appears to give a countable 
infinity of one-parameter families of reflection factors, but this is not so: the 
blocks (x) in (5.5) and (5.6) depend on the boundary parameters only through 
the combination 


2ii9 

TT 


2r] 


TT 


+ 2{r-d) + 1 


(5.11) 


Since {x + 2h) = (x), the reflection factors for d are therefore the same as those 
for d + h. In addition, the freedom to redefine the parameter r] gives an extra 
invariance of the one-parameter families under d ^ 2r — d. Thus the full set of 
options is realised by 

d = 0,l,---r. (5.12) 


Note also that i?^“^(d), the coupling-dependent part of the reflection factor, is 
trivial if d = 0. Thus the limit (5.8) corresponds to r one-parameter families of 
breather reflection factors, and one ‘isolated’ case. This matches the counting of 
conformal boundary conditions (the set of bulk Virasoro primary fields [36]) for 
the Al2,2r+3 minimal models, and the fact that of these boundary conditions, all 
but one have the relevant (^13 boundary operator in their spectra. 

To see how this fits in with our earlier ideas, we swap rj for C, defined by 


2 - 1 ] 

TT 



(5.13) 


Then, starting from (5.5) and relabelling, 


nd) — TT {—h + d + l + C) {—d + l + C) 

^ [d 1 C) [h — d 1 C) 

step 2 
a—1 

= Yl {-d-l-C){-d-l + C){-h + d + l + C){-h + d + l-C) 

1=1—a 
step 2 

a+d—1 

n {-l-C){-l + C){-h + l + C){-h + l-C) . (5.14) 

l=l—a-\-d 
step 2 

Since the first a — d terms in the last product cancel if a > d, it is easily seen that 
coincides with as expressed by (4.10), for the S-matrix (2.9). 
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6 Some aspects of the off-critical ^r-functions 


Our main tool in linking reflection factors to specific boundary conditions will be 
the formula for an exact off-critical gf-function introduced in [11], In this section 
we begin by reviewing the main formula, and then discuss some further properties 
and modifications which will be relevant later. 


6.1 The exact gf-function for diagonal scattering theories 

The ground-state degeneracy, or g-function, was introduced by Affleck and Lud¬ 
wig as a useful characterisation of the boundary conditions in critical quantum 
field theories [8]. In a massive model, the definition can proceed along similar 
lines [9, 2, 11], One starts with the partition function for the the¬ 

ory on a finite cylinder of circumference L, length R and boundary conditions of 
type a at both ends. In the ii-channel description, time runs along the length 
of the cylinder, and the partition function is represented as a sum over an eigen- 
basis {I'lpk)} of the Hamiltonian which propagates states along the 

cylinder: 


^a\d) [^1 (Q;| e 




E 


.-REfRM^L) 


fc =0 


Here I = ML, M is the mass of the lightest particle in the theory, and 

— (Qlbfc) 


( 6 . 1 ) 


( 6 . 2 ) 


where ja) is the (massive) boundary state [26] for the a boundary condition. 
At finite values of I any possible infinite-volume vacuum degeneracy is lifted by 
tunneling effects, making the ground state jV’o) unique. This gives ^(a|o) the 
following leading and next-to-leading behaviour in the large-i? limit; 

lnZ(,|„)[L,f?] ~ + (6.3) 

It is the second, sub-leading term which characterises the massive g'-function. 
Subtracting a linearly-growing piece —the g(-function for the boundary 
condition a at system size I is dehned to be 


lngf |„)(0 = ln ^|°|(0 f\a)L . 


(6.4) 


The constant is equal to the constant (boundary) contribution to the ground- 
state energy Llg*”'’(i?) of the L-channel Hamiltonian iL®*’'‘P(i?) = (R), which 

propagates states living on a segment of length R and boundary conditions a at 
both ends. 

Although this definition is in principle unambiguous, the difficulty of charac¬ 
terising the finite-size states |a) and j'i/’o) limits its utility in practice. An alter¬ 
native expression for g can be obtained by comparing (6.3) with the L-channel 
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representation, a sum over the full set of eigenvalues of ; 

OO 

[L, R] = f2 . (6.5) 

fc =0 

As R is sent to infinity, we have 

ln5|«)(0 = ^ lim L (f2 + 2f\^}L + RE^iM, L)] . ( 6 . 6 ) 

Z it—>oo \ 

I \k=Q / J 

In theories with only massive excitations in the bulk and no infinite-volume 
vacuum degeneracy, lngr|^^(^) tends exponentially to zero at large Z, while in the 
UV limit / ^ 0 it reproduces the value of the sf-function in the unperturbed 
boundary conformal field theory. 

The subleading nature of the term being extracted makes it hard to evaluate 
( 6 . 6 ) directly. In [11], considerations of the infrared (large-Z) asymptotics and a 
conjectured resummation led to the proposal of a general formula for the exact 
^(-function of a purely elastic scattering theory. However in section 6.2 below we 
shall argue that the result of [11] needs to be modified by the introduction of a 
simple symmetry factor whenever there is coexistence of vacua at infinite volume. 
This happens in the low-temperature phases of the Eq, Ej, A and D models. The 
more-general result is 

1 r 

ln5|„>(0 = Inq.) + 7 E / [cp\:'>{e)-6{e)-2<Paam) In (l + + S(Z) 

(6.7) 

where C\a) is the symmetry factor to be discussed shortly, and the functions 
(l)a\o) and (j)ab{G) are related to the bulk and boundary scattering amplitudes 
Sabi9) and R^a\9) by 

= -^^lnR\:\9), ( 6 . 8 ) 

4>ab{9) = . (6.9) 

The boundary-condition-independent piece S(Z) can be expressed in terms of 
solutions to the bulk thermodynamic Bethe ansatz (TBA) equations as 


m = 


iv V - /_ — 


2 ^ ^ n 1 + 1 + 

n=l ai...an=l 

('/’aia2 (^1 T 92)4>a2a3 (^2 ~ ^s) ...</* a„a„+i(^n ^n+1 )) (6.10) 


with 9n+i = 9i, On+i = oi. The functions £a{9)-, called pseudoenergies, solve the 
bulk TBA equations 


£a{9) = MgL cosh 6 > - E / d9'(t)ab{9 - 9')Li,{9') , a = l,...,r 0 , (6.11) 

b=i 
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where La{0) = ln(l + e We also recall that the ground-state energy of the 

theory on a circle of circumference L is given in terms of these pseudoenergies by 

(M, L) =-V / —MaCosheLa{e)+SLM^ (6.12) 

27r 

where SLM"^ is the bulk contribution to the energy. 

In [11] the formula (6.7) was checked in detail, but only for the = 1, = 1 

case corresponding to the Lee-Yang model. Our results below will confirm that 
it holds in more general cases, provided C\a) is appropriately chosen. 

6.2 Models with internal symmetries 

The scaling Lee-Yang model, on which the analysis of [11] was mostly concen¬ 
trated, is a massive integrable quantum field theory with fully diagonal scattering 
and a single vacuum. However many models lack one or both these properties. In 
this section we shall extend our analysis to models possessing, in infinite volume, 
N equivalent vacua but still described by purely elastic scattering theories ^ . For 
the ADE-related theories the N equivalent vacua are related by a global symme¬ 
try Z. We shall see that when these systems are in their low-temperature phases, 
the formula for g{l) as originally proposed in [11] should be slightly corrected, by 
including a ‘symmetry factor’ C\ct)- 

Low-temperature phases are common features of two dimensional magnetic 
spin systems below their critical temperatures. A discrete symmetry Z of the 
Hamiltonian H is spontaneously broken, and a unique ground state is singled 
out from a multiplet of equivalent degenerate vacua. This contrasts with the 
high-temperature phase, where the ground state is ^-invariant. 

The prototype of two-dimensional spin systems is the Ising model, and we 
shall treat this case first. The Ising Hamiltonian for zero external magnetic 
field is invariant under a global spin reversal transformation {Z = Z 2 ), and at 
low temperatures T < Tc this symmetry is spontaneously broken and there is a 
doublet {]+), |—)} of vacuum states transforming into each other under a global 
spin-flip: 

|±)^|T). (6.13) 

The Z 2 symmetry of H is preserved under renormalisation and it also charac¬ 
terizes the continuum field theory version of the model: in the low-temperature 
phase the bulk field theory has two degenerate vacua {]+), |—)} with excitations, 
the massive kinks iL[_|_ ](0), iLj_|-](6*), corresponding to field configurations inter¬ 

polating between these vacua. 

In infinite volume the ground state is either j-|-) or j— ), and the transition from 
j-|-) to j— ) (or vice-versa) can only happen in an infinite interval of time. On the 
other hand in a finite volume V, tunneling is allowed: a kink with finite speed 

Hntegrable models with non-equivalent vacua are typically associated to non-diagonal scat¬ 
tering like, for example, the ^13 perturbations of the minimal A4p.q models [41]. 
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can span the whole volume segment in a finite time interval. We are interested 
in the scaling limit of the Ising model on a finite cylinder and we interpret the 
open-segment direction of length R as the space coordinate. Time is then periodic 
with period L. If the boundary conditions are taken to be fixed, of type at 
both ends of the segment, then the only multi-particle states which can propagate 
have an even number of kinks, of the form 


... Ki_^]{en) , (n even). (6.14) 

For R large the rapidities {6j} are quantised according to the Bethe ansatz equa¬ 
tions 

rsinh6j — ilnR^~^\6j)=7rj, {r = MR, j = 1,2,...) , (6.15) 

where (9) is the amplitude describing the scattering of particles off a wall 
with fixed boundary conditions of type -|- . 

At low temperatures [L, i?] therefore receives contributions only from 

states with an even number of particles. It is conveniently written in the form 
(see for example [42]): 

Z(+|+) [L, R] = iz®, [L, R] + [L, R] (6.16) 

where 

^(+|+)[^, R] = n (^ + (-1)^6-'“*^^'^^) , {I = ML ), (6.17) 

j>0 

i?) is the ground-state energy, and the set {9j} is quantised by (6.15). 
In order to extract the subleading contributions to Z^_|_|_|_^ , as in [11] we first set 

Z(+|+) [L, R] = , (6.18) 

and in the limit R ^ oo use Newton’s approximation to transform sums into 
integrals. The result is 

ln^J5|+)[i,i?] - \ j d0(^cosh(0) + </.l+>(0)-5(0))ln(l + (-l)V'“^^'') 

= 2 In^i®^ (0 - REf^^‘^'\M, L) ,{b = 0,1) (6.19) 

From (6.19) we see that at finite values of / as i? ^ oo 

In [L, R] - In [L, R] ^ +oo (6.20) 

and therefore 

= i (9(S(i))7-“« . (6.21) 
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where L) is the ground state energy for the system 

with periodic boundary conditions. 

Notice that [L, R] takes contributions from states with both even and 

odd number of particles. In fact, following [11], we have 

<7|+>(0 = <7|->(0 = • (6-22) 

In conclusion, the appearance of the extra symmetry factor ~ related 
to the kink selection rule restricting the possible multi-particle states. In the 
Ising model this was seen in the exact formula for the low-temperature partition 
function ^(+|+) [L, R], by its being written as an averaged sum of [R, R] and 

In more general theories with discrete symmetries, similar considerations ap¬ 
ply. Consider a set of i = 1,2,... ‘fixed’ boundary conditions, matching the 
i = 1,2,... N vacua. A derivation of as in [11], including all multiparticle 
states satisfying the Bethe momentum quantisation conditions, will give an incor¬ 
rect answer, since some states will be forbidden by the kink structure. Instead, 
let us construct a new boundary state 


N 

ic') = Ew- 


i=i 


(6.23) 


and consider the partition function Since boundary scattering does not 

mix vacua it is clear that the only effect of replacing |i) by |?7) is to eliminate the 
kink condition on allowed multiparticle states. Thus the counting of states on an 
interval with U at one end and i at the other is exactly the same as it would be 
in a high-temperature phase with the reflection factor at both ends being i?|j^(0), 
and the derivation in [11] goes through to find at large R 

In Zi^uii) [L, R] ~ -RE^^^M, L) + 2 In (/) (6.24) 


where is - up to a linear term - given by (6.7) with (/>|(^^(0) = (j)a\o) and 
C\ci) = 1- On the other hand, since all vacua are related by the discrete symme¬ 
try Z and the finite-volume vacuum state |'0o) niust be symmetrical under this 
symmetry, (V'olO = (V’ojj) Vi, j, and so 

Hence 

5|i>(0 = ^5(0 (6-26) 

and g\i){l) is given by the formula (6.7) with C\i-^ = 
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An immediate consequence is that in 
Instead, 

lim = 

/—>00 


the infrared, g\i){l) does not tend to one. 




1 

y/N ■ 


(6.27) 


This perhaps-surprising claim can be justified independently. In infinite volume, 
boundary states can be described using a basis of scattering states [26]: 


|B> = (i +1 /“ K"’’(e)Aj-e)Am +... j |o) 


(6.28) 


where |0) is the bulk vacuum for an infinite line. In finite but large volumes, the 
same expression provides a good approximation to the boundary state, the only 
modification normally being that the momenta of the multiparticle states must 
be quantised by the relevant Bethe ansatz equations. However a subtlety arises 
when there is a degeneracy among the bulk vacua, in situations where the bound¬ 
ary condition distinguishes between these vacua. (Similar issues are encountered 
in comparing exact bulk VEVs of fields with those obtained from finite-volume 
approximations [43, 44].) Take the case of a fixed boundary condition which picks 
out one of the degenerate vacua, say i: then the state jO) on the RHS of (6.28) 
is ji). However when calculating the finite-volume ^(-function, we must consider 
where {'ipol is the finite-volume vacuum. (Note that the bulk energy is 
normalised to zero when considering (6.28), so this inner product will give us g 
directly, rather than Q.) In large but finite volumes, the tunneling amplitude 
between the infinite-volume bulk vacua is non-zero and so the appropriately- 
normalised finite volume ground state is the symmetric combination 


(V'ol 


N 


.7 = 1 




(6.29) 


The limiting value of g in the infrared is therefore not 1, but 1/^/N, as found in 
the exact calculation earlier. 

Table 1 lists the central charge, the Coxeter number, the symmetry group Z, 
and the number of degenerate vacua for the ADET models. For the Ar, Dr, Eq 
and El theories at low temperatures there is a coexistence of N vacuum states 
ji), i = 1,2,..., N, with N = r 1, 4, 3 and 2 respectively. The corresponding 
symmetry factor C\i'^ , for fixed-type boundary conditions which single out a single 
bulk vacuum, is always 1/^/N. (For more general boundary conditions one can 
anticipate other symmetry factors, but we shall leave a detailed discussion of this 
point for another time.) 


6.3 Further properties of the exact gf-function 

An important property of the sets of TBA equations that we are considering 
concerns the so-called Y-functions [37], 

Yj(0) = (6.30) 
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Model 

Ceff (0) 

Coxeter number h 

Global symmetry Z 

N 


2r 

r+3 

r + 1 

Zr_l_i 

r -b 1 

Dr {r even) 

1 

2r-2 

Z 2 X Z 2 

4 

Dr {r odd) 

1 

2r-2 

Z 4 

4 

Eg 

6 

7 

12 

Z 3 

3 

Et 

7 

10 

18 

Z,2 

2 

Eg 

1 

2 

30 

Wji 

1 

Tr 

2r 

2r-|-3 

2r -|- 1 


1 


Table 1; Data for the ADET purely elastic scattering theories. The symmetry factor 
C\i) for fixed-type boundary conditions is 1/^fN in each case, where N is the number of 
degenerate vacua in the low-temperature phase. 


These satisfy a set of functional relations called a Y-system: 

- ^) = n (1 + , (6.31) 

6=1 

where is the incidence matrix of the Dynkin diagram G labelling the TEA 
system. Defining an associated set of T-functions through the relations 

[G] 

Ya{e) = , i + Ya{e) = Ta{e + ii)Ta{e-ii) (6.32) 

c=l 


we have 

[G] 

Ta {9 + ii)Ta{e - if) n (re(0))-^“ = 1 + Y-\e ). 

c=l 

In addition the T-functions satisfy 


Ta{e + - !^) = 1 + J] [Tmr 


ZTT, 

T 


6=1 


(6.33) 


(6.34) 
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Fourier transforming the logarithm of equation (6.33), solving taking the large 6 
asymptotic into account, and transforming back recovers the (standard) formula 


T M ^ r 

\nTa{9) = - - ^ cosh ed9'xab{0-0')Lb{O'), a = l,...rg, (6.35) 

^ cos ^ Jk 

where [37, 17] 

Xafe(0) =(2cosh(/c7r//i)l - =-^^5fb(0). (6.36) 

The result (6.35) allows us to make a connection between the exact ^(-functions 
for our parameter-dependent reflection factors and the T functions. Taking the 
reflection factors defined by (4.8) and (4.6) 

=R^^\9)/{S^b{0-ilC)S^b{0 + ilC)) (6.37) 

and using (6.7), we find 

1 r 

lng^b,c)il) = d9'xabi9'-ilC)Lbi9') 

"bE [ d9'xab{9' + ilC)Lb{9'). (6.38) 

^ a=l dK 

Comparing this result with (6.35) and using the property Tb{9) = Tb{—9), 

In5-16,0(0 = ln5'(0 + Inrb(f^C') - cos(^C') (6.39) 

or, subtracting the linear contribution; 

^|?o(0 = ^^°^(0r6(i^C'). (6.40) 

Exact relations between g- and T- functions in various situations where the bulk 
remains critical were observed in [38, 39, 40]. These were extended off-criticality 
to a relation between the Q- and T- function of the Lee-Yang model in [10]. 
The generalisation of [10] proposed here relies on the specific forms of our one- 
parameter families of reflection factors, and thus provides some further motivation 
for their introduction. 

It is sometimes helpful to have an alternative representation for the infinite 
sum S(f) in (6.7), the piece of the ^-function which does not depend on the specific 
boundary condition. It is readily checked that, for values of I such that the sum 
converges, 

OO ^ OO - 

= E - E 

n=l n=l 
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where 


Kabi0,6') = - 


and 


The identity 


2 

Hab{0,e') = 


{^ab{G + o') + (ftahiP — O'^ 


y/l + e6i,(d') 


-4ab{0 - o' 


\/1 + \/1 + 


(6.42) 


(6.43) 


- 

lnDet(/ - M) = Trln(/ - M) = - V -TrM*" 

^ n 

n=l 


(6.44) 


then allows S(^) to be rewritten as 

S(0 = i lnDet(^^—. (6.45) 

This formula makes sense even when the original sum diverges, a fact that will 
be used in the next section. 


7 UV values of the ^r-function 

Taking into account the effect of vacuum degeneracies described above, a first 
test of the reflection factors found earlier is to calculate the I 0 limit of (6.7). 
The value of 5 '|q,)( 0) should match the value of a conformal field theory g'-function, 
either of a Cardy state or of a superposition of such states. 

As / ^ 0, the pseudoenergies £a{0) tend to constants, which we shall denote as 
Cq. Their values were tabulated by Klassen and Melzer in [7], who also observed 
an elegant formula the integrals of the logarithmic derivatives of the bulk S-matrix 
elements, later proved in [15] : 


d6 4>ab{0) = -Mab , 


(7.1) 


where Nab is related to the Cartan matrix C of the associated Lie algebra by 
M = 2C“^ — 1.^ We shall also need the integrals of the logarithmic derivatives 
of the reflection factors. Writing 


-Ra= n 

x&A 


( ' 

-1 2 -1 


^The Cartan matrix for Tr is taken to be 




I i j 


(7.2) 
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for some set A, the required integrals are 



-2 ^ - 0 sign[x//i] 

xeA 


(7.3) 


with sign[0] = 0. For the minimal ADE reflection factors found in section 3 these 
evaluate to 

f d0cf>l^>(0) = l-Afaa, (7.4) 

jR 

while for the reflection factors (5.4), 



(7.5) 


To calculate the UV limit of S(/), we dehne a matrix A4 whose elements are 


■Mab = 


Af2b 

1 + 6*^“ ’ 


(7.6) 


in terms of which 


^(0) = ) E -« + ■■■+'”) = -I -«)■■■ (1 - 

4 n 4 

n=l 


where ei,..., are the eigenvalues of A4. When one of these eigenvalues is larger 
than 1, the infinite sum does not converge (this is the case for Ar, r > 5, D^, 
r > 4 and E^.) but the RHS of (7.7) gives the correct analytic continuation, since 
it follows from (6.45). 


7.1 Tr 


We start with the theories, whose UV limits are the non-unitary minimal 
models A42,2r+3- There are r + 1 ‘pure’ conformal boundary conditions, and the 
corresponding values of the conformal ^r-functions can be found using the formula 


9{l,l+d.) 


*S'(l,l+r);(l,l+d) 
l*S'(l,l+r),(1,1)1^/^ 


d = 0,..., r 


(7.8) 


where (1,1 + r) is the ground state of the bulk theory, with lowest conformal 
weight, (1,1) is the ‘conformal vacuum’ ~ with conformal weight 0 - and S is 
the modular S-matrix (see for example [45]). The components of S for a general 
minimal model Mp'p are 


S{n,my,ip,a) = sin(7r^np) sin(7r—mcj) , (7.9) 

where 1 < n, p < p' — 1, 1 < m,a < p — 1 and, to avoid double-counting, m < 
and p < For the Tr theories, p' = 2 and p = 2r -|- 3, and so n and p are both 
equal to 1 while m and a range from 1 to r -|- 1, as in (7.8). 
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The values predicted by (7.8) should be compared with the values of 5 '|q)( 0) 
calculated from (6.7). For the boundary-parameter-independent part of the re¬ 
flection factor from (5.4), which is the d = 0 case of the ‘quantum group reduced’ 
options (5.12), the second term of (6.7) evaluates to 


1 ^ r 

4 E / - m - 2(t>aa{2e)) In (1 + e-^“) 

a=l 

j . / _:._2 / TT \ 

= (l + —T 


2r+3 


= In 


while S(0), calculated using (7.7), is 



--ln((l-ei)...(l-e,.)) = ^ 


-In f 


V\/2f+3 


sm 


(r-|-l)7r 

2r-|-3 


(7.10) 


(7.11) 


Adding these terms together and using some simple trigonometric identities re¬ 
veals a dramatic simplification: 


for all r. This suggests that the minimal Tj. reflection factors (5.5) describe bulk 
perturbations of the boundary conformal field theory with (1,1) boundary con¬ 
ditions. Notice that among all of the possible conformal boundary conditions 
in the unperturbed theory, this is the only one with no relevant boundary op¬ 
erators, matching the fact that the minimal reflection factors (5.5) have no free 
parameters. 


7.2 Aj., Dr and Er 

The ADE cases exhibit an interestingly uniform structure; substituting (7.1) and 
(7.4) into (6.7) shows that the UV limit of the second, reflection-factor-dependent, 
term of (6.7) is always zero when the minimal reflection factors from section 3 are 
used. This result can be confirmed by examining the contributions of the blocks 
(x) of S{29), and (x) (or (x)) of i?l"^(6*), as follows. Since the bulk S-matrix can 
be written as 

Sabie) = n {^}(^) (7-13) 

where {x} = (x -|- 1) (x — 1), Sab{29) can be written similarly; 

Sab{29) = n N(^) (7-14) 

X^Aab 
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where 


[x] 


— 1^ ^2/i — X + —X —^ 


Summing the contributions of each block [x], it is easily seen that 

'2x 

ceAoa ^ 


2 / decPaam= V 

J M A 


h 


+ Sx.i — 2 


(7.15) 


(7.16) 


The minimal reflection factors discussed in section 3 can be written in a similar 
way; 

Ra= W fx (7.17) 

X&Aaa 

where each fx is either (x) or (x). The contribution to f d0 from each (x) 

is —2 + 2x//i + 2(5x,i and from each (x) is —2 + 2x/h. Noting that every minimal 
reflection factor contains the block (1) exactly once, we find 

dect>\^\9)= Y. (^^ + 24,1-2), (7.18) 



and so 

jje ( 4 “>( 0 ) - 5{e) - 2ci>aam) = 0 


(7.19) 


for every A, D and E theory, as claimed. Working backwards, this gives a general 
proof of the formula (7.4), given (7.1). 

With the reflection-factor-dependent term giving zero, the sum of S(0) and 
the symmetry factor should correspond to a conformal (/-function value. To find 
these values in a uniform way, we shall use the diagonal coset description gi xgi/02 
where Qi is the affine Lie algebra at level I associated to one of the A, D or E 
Lie algebras [46]. Coset fields are specified by triples of g weights at 

levels 1, 1 and 2, and the (/-function for the corresponding conformal boundary 
condition can be written in terms of the modular S-matrix of the coset model 




'S'{0,0;0}{0,0;0} 


(7.20) 


The integrable highest weights p of g can be expressed in terms of the fundamental 
weights Aj of g as: 

re 

p = Y'^i-^i ('^• 21 ) 

i=0 


where the non-negative integers n*, i = 0, ...,rg, are the Dynkin labels. The 
fundamental weights can be written as 


Ao = (0,1,0) (7.22) 

0 ) 
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where A* are the fundamental weights of g and the colabels, Oj, are given for each 
case below: 


Ar :di = 1 for all i 

Dr : di = dr-i = dr = 1 all others are 2 

^6 : a* = (1,1, 2, 2, 2,3) (7.23) 

Ey ;ai = (1,2,2,2,3,3,4) 

Eii : di = (2, 2, 3,3,4,4, 5,6). 


The highest weights /ig of the corresponding finite-dimensional Lie algebra, g, 
can be written in a similar way, with Dynkin labels, ni,..., rir ^, taken from the 
same set: 

^■b 

t^o = '^niXi. (7.24) 

i=l 

Once the level I of g has been specified, the Dynkin labels must satisfy 

^•b 

I = no + nidi (7.25) 

i=l 


so the representation of Qi is completely determined by the Dynkin labels of the 
corresponding representation of g. For example, the label 0 is given to the vacuum 
representation, where n* = 0, i = 1,..., for both gi and Q 2 - 

In appendix A we show among other things that the coset conformal g- 
function (7.20) depends only on the level 2 weight, via the level 2 modular S- 
matrixll: 


o(2) 

_ _ 
-gp- 

V *^00 


(7.26) 


Clearly, this will not fix the coset field in general. For example, the number of 
boundary conditions with ^r-function equal to go is r -|- 1, 4, 3, 2, 2 for the Ar, 
Dr, Eq, Ei and Eg models respectively. More details on this point are given in 
appendix A. 

Algorithms for computing the modular S-matrices are given by Gannon in [47]; 
Schellekens has also produced a useful program for their calculation [48]. The level 
2 modular S-matrix elements for A, D and E theories needed to calculate the UV 
values of the g-functions using eq.(7.26) are given in table 2. The representations 
are labelled by the Dynkin labels, n*, i = 0,... ,r. The number of coset fields 
corresponding to each label is equal to the order of the orbit of that level 2 weight, 
under the outer automorphism group 0(g). Note for Dr, r even, the weights pr-i 
and pr are in different orbits, each with order 2, whereas for r odd they are in the 
same orbit with order 4, so in both cases there are 4 coset fields with the same 
g-function value. 


I^We’d like to thank Daniel Roggenkamp for his help in explaining this observation. 
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S-matrix element 

Labels 


~ ih+ 2 )Vh ® {h+ 2 ) llfc=l Sin ) 

c(2) _ 2^ -in r '"in \ 

*OPj - (?,+2)vdi ® V '*+2 

0 : no = 2 , rii = 0 for i = 1,..., r 

1 f 0 * J 

p - no = 1, ni = 1 ^ ^ 

for j = 1,..., [h/ 2 ], 

where [x] is the integer part of x 
and h is the Coxeter number. 

Dr 

o(2) 1 h 

^00 — A\J r 

0(2) _ c(2) _ oc(2) 

q{2) _ 0(2) _ 1 

— 2^2 

0 : no = 2 , m = 0 for i = 1,..., r 

Pi : no = ni = 1 , all other m = 0 

A f 0 * J 

p - no = 0 , ni = 1 ^ i = J 

for j = 2,...,r/2 

Pr--! : no = nr-i = 1, all other ni = 0 
Pr ■ no = nr = 1, all other ni = 0 

Ee 

^00 ^ ^ (i^) 

= ^sin( ) forj = l,2 

0 : no = 2 , ni = 0 for t = 1 , 1 , 2 , 3 , 3,4 

Pi : no = 1, ni = 1 all other ni = 0 

P2 : no = 0, n2 = 1 all other Ui = 0 

Er 

c(2) _ 2 ^in f 47 r \ 

*^00 “ v^ 7 i+ 2 ®^’^ V^-r2^ 
<\=2^sin(^) 

S^ol = J+ 2 ^^^{h%) 

c(2) _ 2,/Z^sin 

h +2 \^h-|-2 J 

0 : no = 2 , ni = 0 for t = 1 ,..., 7 

Pi : no = 1 , ni = 1 all other m = Q 

p2 : no = 0, n2 = 1 all other Ui = 0 

P3 : no = 0, no = 1 all other m = Q 

Es 

c(2) _ 0(2) _ 1 

*^00 “ ^0P2 ~ 2 

c{2) _ 1 

•^OPi “ s /2 

0 : no = 2 , ni = 0 , for all i = 1,..., 8 
P2 : no = 0, n2 = 1, all other ni = 0 

Pi : no = 0 , ni = 1 all other ni = 0 


Table 2: Level 2 modular S-matrix elements for A, D and E models 
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For each coset, the possible CFT values of the g'-function can now be calculated 
using eq.(7.26) and the modular S-matrix elements given in table 2. Working case- 
by-case, we checked these numbers against the sums h\C\a) +^(0), the UV limits 
of the off-critical gf-functions g{l) for the minimal reflection factors described in 
section 3. In every case, we found that 

lng{l)\i=o = InCiQ) -h i;(0) = In 5-0 (7.27) 

provided that the symmetry factors C\a) were assigned as in table 1. The explicit 
^(-function values are given in table 3. 



C\a) 

90 

Number of fields 


1 

Vr+l 

( ih+2)Vh {hl2) nU ^ sm 

r -|- 1 

Dr 

1 

2 

1 /2U/4 

2 \r) 

4 

Eg 

1 

U 3 


3 

Ej 

1 

U2 


2 

00 

1 

i 

_VI_ 

2 


Table 3; UV ^-function values calculated from the minimal reflection factors for the A, 
D and E models 


For some minimal models, we can compare go to the values of the < 7 -function 
corresponding to known cases, thereby matching our reflection factors to physical 
boundary conditions. For the three-state Potts model {A 2 ), the tricritical Ising 
model (Ej) and the Ising model {Eg) the corresponding boundary condition is 
the ‘fixed’ condition in each case [49, 50, 8 ]. Note the number of coset fields with 
^(-function equal to 5^0 corresponds to the number of degenerate vacua, and hence 
to the number of possible ‘fixed’ boundary conditions for all A, D and E models. 

8 Checks in conformal perturbation theory 

The off-critical (/-functions only match boundary conformal field theory values 
in the far ultraviolet. Moving away from this point one expects a variety of 
corrections, some of which were analysed using conformal perturbation theory in 
[10]. The expansion provided by our exact (/-function result is instead about the 
infrared, but convergence is sufficiently fast that the first few terms of the UV 
expansion can be extracted numerically, allowing a comparison with conformal 
perturbation theory to be made. In [11] this was done for the boundary Lee- 
Yang model; in this section we test our more general proposals for the case of the 
three-state Potts model. Since the treatment of conformal perturbation theory 
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in [10] concentrated on the non-unitary Lee-Yang case, we begin with a general 
discussion of the leading bulk-induced correction to the g'-function in the (simpler 
to treat) unitary cases. 

Consider a unitary conformal field theory on a circle of circumference L, per- 


turbed by a bulk spinless primary field (p with scaling dimension 

The perturbed Hamiltonian is then 

— -|- A(^ . 


H = Ho + A#i 

(8.1) 

where 


(8.2) 

and 


(8.3) 


For A real in the ADE models we also expect a A-M relation of the form 


M(A) = ^t|A|i/(2-x^) ; |A(M)| = f--j (8.4) 

with K a model-dependent constant. 

Leaving the boundary unperturbed, we have a conformal boundary condition 
a, with boundary state |a). Set g\a) = = (a|0) and = {a\^p), where |0) 

and \ip) are the states corresponding to the fields 1 and (p. Since the theory is 
unitary, |0) is also the unperturbed ground state; and since p is primary, we have 
( 0 |(^| 0 ) = 0 . 

The aim is to calculate ln^|Q,^(A,L) = ln(a|n) where |a) is the unperturbed 
CFT boundary state, and jfl) is the PCFT vacuum. This will be a power series 
in the dimensionless quantity ; here, we shall obtain the coefficient of the 

linear term, The calculation follows [10], but is a little different (and in fact 
simpler) because the theory is unitary, so that the ground state is the conformal 
vacuum |0). 

First-order perturbation theory implies 

|A1) = |0) -|- A flajV'a) + • • • (8-5) 

a 


where the sum is over all states excluding |0) , (V'ajO) = 0 and 


(0|#o|0) - (V^al^olV’a) ■ 


( 8 . 6 ) 


Since the theory is unitary, (0|iLo|0) - (V'aj^^ojV'a) = -(V'alx(-^o + -^^o)|V’a)- Using 
rotational invariance as well. 


f^alV’a) = 


AL^ W |V>a)(V;a|V9(l)|0) 

(2vr)l-^ 4" (V’ajUo+IolV'a) ■ 


(8.7) 


28 



Hence 


XL 


‘l — Xu 


1 


|Si> = |0)-p^(l-P);^(l-FMl)|0)+.,, 


m 

where P = |0)(0| is the projector onto the ground state. Using the formula 

\T2-Xtp l-l J _ 

(a|U) = (a|0) - (l-P) -j (?^°+^°(l-P)^(l)|0) + ... (8.9) 


Lq+Lq 


Since (0|(/9|0) = 0 and ^°|0) = g*‘^|0) this last 

expression simplifies to 

\t2—Xip rl 

(«|U) = (a|0) - / dqq^‘^-^{a\ip{q)\0) + ... (8.10) 

{ 2 TTy Jq 

Now (a|(^(g)|0) is a disc amplitude, and by Mobius invariance it is given by 

{a\^{q)\0) = {1 - . (8.11) 

(This is a significant simplification over the nonunitary case discussed in [10], 
where the corresponding amplitude had to be expressed in terms of hypergeo¬ 
metric functions.) 

Taking logarithms, 

\ t2-x^ gf j-l 

ln^|Q)(A,L) = lng|„)-—— - — dq q^‘^~^{l - + ... 

(2vr)^ g\a) Jo 

= ln5|„)+4“^L2-a;^ + ... (8.12) 


and doing the integral, 


d 


l“> 

1 


1 g|«) 

2(27r)i-^^ 5|a> 


H (1 Xip^Xipl^'j 


(8.13) 


where B{x,y) = r(x)r(y)/r(x-|-y) is the Euler beta function. 

This is a general result. As a non-trivial check we specialise to the 3-state Potts 
model, described by the A 2 scattering theory. There are three possible values A, B 
and C of the microscopic spin variable, related by an S 3 symmetry. At criticality 
the model corresponds to a c = 4/5 conformal field theory. The primary fields 
are the identity I, a doublet of fields {'0,V’^} of dimensions = 2/3, the 

energy operator e of dimensions A^ = A^ = 2/5 and a second doublet of fields 
{a, } with dimensions A^- = A^- = 1/15. The bulk perturbing operator ip which 

leads to the A 2 scattering theory is e, and so x^ = = 4/5. Boundary 

conditions and states for the unperturbed model are discussed in [55, 49]. One of 
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the three ‘fixed’ boundary states, say |^), can be written in terms of Vha-Ishibashi 
states as [55] 

\A) = \i) = K[\I)) + X\E)) + \i;)) + \^l;^)) + X\a)) + X\a^))] (8.14) 

where = (5 — \/5)/30 and = (1 + \/5)/2. Hence 

In = In X =-0.5961357674... , = X = 1.2720196495 ... (8.15) 

Putting everything into (8.13), the CPT prediction for the coefficient of the first 
perturbative correction to the ^-function for fixed boundary conditions in the 
three-state Potts model is 

= -3.011357884... (8.16) 

The 3-state Potts model also admits ‘mixed’ boundary conditions AB, BC and 
CA [56]. For later use we recall from [55] that the corresponding boundary state 
is 

\AB) = K [xV)) - ^“V)) + )) - ^“V)) - ^“^1 )) 

(8.17) 

so that, for the AB boundary, we find instead 

= (8.18) 

The results (8.15) and (8.16) can be compared to the numerical evaluation of 
our exact ^f-function result, written in terms of AT®/® using Fateev’s formula [51] 

« = ^J|^(27r)^/®(7(2/5)7(4/5))'^/i2 ^ 4.504307868 ... (8.19) 

where 7(x) = F(x)/F(l — x). Setting x = = —{1/ ^)6/5 (i-eixiembering that 

(T — Tc) oc A < 0) a fit to our numerical data yielded 

In g{l) = -0.596135768 - 0.49999991 1 - 3.0113570 x 

- 0.909937x2 0.3982 x^ -M.Ox"^ ... (8.20) 

which agrees well with (8.15) and (8.16). The match of the constant term in 
(8.20) with Ing'i^i^ from (8.15) is guaranteed by our exact formula, and so serves 
as a check on the accuracy of our numerics. A calculation of the coefficient of the 
irregular (in x) term, proportional to I, as in [10], predicts the value 0.5, again in 
good agreement with our numerical results. 

9 One-parameter families and RG flows 

Just as the minimal reflection factors were tested in section 7, we do the same 
now with the one-parameter families of reflection factors. These reflection factors, 
R^a’^\ depend on the parameter C and are given in (4.6) and (4.10). 
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9.1 The ultraviolet limit 

If we use the reflection factors as input to calculate the ^-function, eq. (6.7), and 
take the limit I 0 then 


9^'^^ ^ 9\d,C){^)\i=o — goTd , (9.1) 

where = Td{9)\i=o is a 0-independent constant. From (6.32) and (6.30) we 
have 

Td = VTT^ (9.2) 

where the values can be found in [7]. For every ADET theory, eq. (9.1) leads 
to a possible CFT g'-function value. 

For M > 0 there will be many massive bulk flows as the parameter C is varied. 
However, it should be possible to tune C, as the limit ^ ^ 0 is taken, so as to give 
a massless boundary flow between the conformal ^(-function corresponding 
to the UV limit of the reflection factor and g, corresponding to the UV 

limit of the minimal reflection factor. These flows are depicted in figure 2. Note 
that the UV sf-function corresponding to the minimal reflection factors {go for 
the ADE cases and 5(1,1) for Tj.) is the smallest conformal value in all cases so, by 
Affleck and Ludwig’s 5-theorem [8], this is a stable fixed point of the boundary 
RG flow. 


Massive Bulk Flow 


Massless 
Boundary 
Flow 



For the R. case we find 

5'"*^ = 9{i,d+i) foi' d = l,...,r. (9.3) 

This is consistent with a simple pattern of flows 

(1,d-|-l) ^ (1,1) for d = l,...,r. (9-4) 
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The conformal values corresponding to the UV limits of for the A, D and 
E cases are given in tables 4 and 5. Again we expect there to be boundary flows 


go (9.5) 

in each case. Notice that in many cases the UV ^r-function values are sums of ‘sim¬ 
ple’ GET values, corresponding to flows from superpositions of Cardy boundary 
conditions, driven by boundary-changing operators. 



Dr 

gW = g\^} = gp^ 

5,|2> = g\r-l) = g^^ 

g,k/2> — g\r/2+i) _ g^^^^ for r even 
= 9ph/2 ^ 

grid = (i + l)gQ for i = 1 ,..., r - 2 

g\r-l) ^ g\r) ^ ^ g^^ 


Table 4: UV g-function values for A and D models 



Eo 

E'j 

Eg 


gW 

gpi 

gpi 

go + dpi 

yPT 


go + gp 2 

go + gp 2 

2 do + dpi 

5 I 3 ) 

gm 

9pi + gp2 

gpi + dps 

2 do + “^gpi 

g\^) 


go + 5pi + “^gpi 

do + ‘^gp2 

3do + ‘^gpi 

yiw 


go + ^gp 2 

5do + 3(7pj 

yisr 


‘^gpi + ‘^gps, 

5do + 

yiTT 


3do + ^9p2 

9do + Qgpi 

yisr 



16go + 12gp^ 


Table 5: UV ^-function values for Eq, Er and Eg models 


The conjectured flow for the three-state Potts model (A2), gp^ go, corre¬ 
sponds to the ‘mixed-to-fixed’ flow (AB A) found by Affleck, Oshikawa and 
Saleur [49], and by Eredenhagen [52]. 

Similarly, the flow gp^ go in the tricritical Ising model (E7) corresponds to 
the ‘degenerate-to-fixed’ ((d) ^ (—) or (d) ^ (-|-)) flows of [53, 52]. Notice that 
go + dpa ^ do in Ey matches the GET g values for the flow (— ) © (0+) ^ (— ) 
conjectured in [52]. However this latter flow is driven by the boundary field with 
scaling dimension 3/5, which is inappropriate for the E7 coset description. A 
more careful analysis shows that our flow, which is driven by a boundary field 
with scaling dimension 1/10, must start from either (+) © (0+) or (—) © (—0), 
and flow to (+) or (—). This serves as a useful reminder that the sf-function 
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values alone do not pin down a boundary condition, and that this ambiguity can 
be physically significant in situations involving superpositions of boundaries. 


9.2 On the relationship between the UV and IR parameters 

The one-parameter families of boundary scattering theories introduced above 
should describe simultaneous perturbations of boundary conformal field theories 
by relevant bulk and boundary operators. The action is 

(9.6) 

where is the unperturbed boundary CFT action. We suppose that the 

boundary condition is imposed at x = 0; in general it might correspond to a 
superposition of n\a) Cardy states. Denoting these by |c), c = 1,2,... ,n\a)-, the 
boundary is in the state 

"■la) 

I") = ’ (9-7) 

C=1 

with n|c) G N. The bulk perturbing term is 

/ O roc 

dx / dy (p{x,y) (9.8) 

-oo J —oo 

while the boundary perturbing part is 

^BND _ ,, 

- 2 ^ 

c^d=l 

The operators 4>(^c\c) hve on a single Cardy boundary, while the 4’(^c\d) with c / d 
are boundary changing operators. Since ip^dd) = ^\d\c)^ ™ ^ unitary theory we 
also expect (see for example [54]) 


dy (l>{c\d){y) ■ 


(9.9) 


y{c\d) ^((i|c) • (9.10) 

Using periodicity arguments to analyse the behaviour of the ground-state 
energy on a strip with perturbed boundaries as in [37] and [2] (see also (9.15) and 
(9.16) below), one can argue that the scaling dimensions of the fields 4>(^c\d) must, 
in these integrable cases, be half that of the bulk perturbing operator ip : 


for the ADE systems and 



2 

h + 2 


Xin 2 — h 


(9.11) 


(9.12) 
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for the Tr models. The results recorded in eq. (9.3) and tables 4 and 5 provide 
information about the conformal boundary conditions associated with the one- 
parameter families of reflection factors. We see that for and the boundary is 
always in a pure Cardy state, while for Dj., Eq and Ef this is true only in one case 
per model, and in the iHs-related theories the UV boundary always corresponds 
to a non-trivial superposition of the states |±) and |free). This observation fits 
nicely with the conformal field theory results for the Ising model [55]: from (9.11) 
we see that for the Eg coset description the integrable boundary perturbation 
must have dimension = 1/16, and indeed the only boundary operators with 
this dimension in the Ising model are </(±|free} and i/>(free|±> • 

For simplicity, we shall only discuss the cases involving a single Cardy bound¬ 
ary, where 

/ OO 

dy (t){y) . (9.13) 

-OO 

As in the Ising and Lee-Yang examples of [26, 2], a simple formula is expected to 
link the couplings A and /x of bulk and boundary fields to the parameter C in the 
reflection factors. However, without a precise identification of the operator cj) it 
is hard to see how such relation can be determined. Even so, a general argument 
combined with a numerically-supported conjecture allows the relation formula to 
be fixed up to a single overall dimensionless constant. This goes as follows. From 
section 6.3 we know that in all cases 

4°C>(0 = ^^°H07^d(*fC') (9.14) 

where is the TBA-related T-function, and is the CPT ^-function corre¬ 

sponding to the minimal reflection factor, for which there is no boundary pertur¬ 
bation. In addition, Td{9) = Td{6,l) is even in 6 , Td{9) = Td{—9), and periodic, 
Td{9 + iTT^^) = Td{9), and so it can be Fourier expanded as 

TdiifC) = Coil)+ ^^Ck{l) cos . (9.15) 

We can now use the observation of [24] that Ta{9,l) admits an expansion with 

2 h 

finite domain of convergence in the pair of variables a± = (/e^®)^+2 to see that 

4 h 2 h 4 h 

Co{l) = Co -|- 0{l ^+2 ) , Clil) = Cl I ^+2 -|- 0(1 ^+2 ) . (9.16) 

We also know that the minimal ^-function has an expansion 

CXD 

In (/) = In 

k=l 

while the conformal perturbation theory expansion of form (see [10]) 

OO 

lng®^^(A,^,L)= c^4^l1-"^)”^(AL2-2-<^)T (9.18) 

m,n=l 


34 



Comparing (9.14) - (9.17) with (9.18) we conclude that, so long as cio / 0, the 
relationship between C and ^ must have the form 

/ Ott \ f 2 tt \ 

^i = lio cos Mh +2 = 'jlQ cos (9.19) 

where /Uq is an unknown dimensionless constant. However the result (9.19) can 
only be correct if 1 — = 2h/{h + 2). This is true only in the non-unitary 

models, and indeed it reproduces the Lee-Yang result of [11] when specialised to 
Ti. For the ADE theories, l — x^ = h/{h+2) and we conclude that cio, the first fi- 
dependent correction to G, must be zero. (This is not surprising since in a unitary 
CFT this correction is proportional to ( (j) = 0.) The first contribution is then 

at order = 0(M^“^*</>), and at this order there is an overlap between the 

expansions of Td{6,l) and This leads to the less-restricted result 


H = koiz + cos 


M ^+2 = k^iz — cos 


h + 2 


(9.20) 


where now both feg and z are unknown constants. If we now consider the Ising 
model, then the n — C formula is known [26]. Written in terms of C it becomes 

= 2M{1 + cos(f C)) ^ /i = 2\/Mcos(f C). (9.21) 

Thus the boundary magnetic field is an even function of C. It is then tempting to 
conjecture that 2" = 1 for all the g\i^c) cases in the Aj. models, and, to preserve the 
perfect square property, that z is either 1 or —1 in all other ADE single boundary 
condition situations; 


1 : = AC--. „ ^ cos (^); 


(9.22) 

(9.23) 


We now return to the physical picture of flows parametrised by C depicted in 
figure 2. At A = 0 the bulk mass is zero, and the only scale in the problem is that 
induced by the boundary coupling g. The massless boundary flow down the left- 
hand edge of the diagram therefore corresponds to varying from 0 to oo. If A 
is instead kept finite and nonzero while J/uj is sent to infinity, the flow will collapse 
onto the lower edge of the diagram, flowing from go in the UV. Eor the ^r-function 
calculations to reproduce this behaviour, the reflection factor should therefore 
reduce to its minimal version as ^ oo. Eor (9.22), \g\ oo corresponds to 
C ioo, which does indeed reduce the reflection factor as required. On the 
other hand, taking j/xj ^ oo in (9.23), requires C {h+ 2)/2 + ioo. While the 
reduction is again achieved in the limit, real analyticity of the reflection factors is 
lost at intermediate values of g. Eor this reason option (9.22) might be favoured, 
but more detailed work will be needed to make this a definitive conclusion. 
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In fact, the proposal (9.22) can be checked at /x = 0 in the 3-state Potts model 
{h = 3), as follows. Consider the results (8.16) and (8.18) and set 

di = ^1 + d^^^ = -0.683763720... (9.24) 

According to the conclusions of section 9.1 and eq. (9.22), 

= af°}(0T’i(fC)|c=5/2 , (9.25) 

and di should match the coefficients ti of in the expansion of the function 
Ti{6,1)\q=q about I = 0. Noticing that Ti{9,l) = T- 2 { 6 , 1 ) = Ti,y{GJ) we can use 
table 6 of [10]: lnrLY(^7r(6 -|- 3)/6) = e{i7r{b + 3)/6), C = 5/2 corresponds to 
6 = 2, = he = —0.6852899839, and we find 

ti ~ 0.9977728224 he = -0.683763721... (9.26) 

which within numerical accuracy is equal to di. 

The study of the cases where the theory is perturbed by boundary changing 
operators is more complicated. Our results are very preliminary, and we have 
decided to postpone the discussion to another occasion. 

10 Conclusions 

Since the initial studies of integrable boundary scattering theories, surprisingly 
little progress has been made on the identification of the many known boundary 
reflection factors with particular perturbations of conformal field theories. One 
aim of this paper has been to use the recently-discovered exact expression for the 
ground-state degeneracy g to begin to fill this gap. Working within the framework 
of the minimal ADET models we were able to identify a special class of reflection 
factors possessing a very simple conformal field theory interpretation: in the 
ultraviolet limit they describe the reflection of massless particles against a wall 
with fixed boundary conditions |z) matching the low-temperature vacua. 

A collection of one-parameter families of amplitudes was also proposed, and 
partially justified by a suggested quantum group reduction of the boundary sine- 
Gordon model at = Ibvr/(2r -|- 3). 

The consistency of these novel families of boundary scattering theories was 
supported by our ^(-function calculations: they appear to correspond to pertur¬ 
bations of Cardy boundary states, or of superpositions of such states, by relevant 
boundary operators, with each boundary flow ending at a fixed boundary condi¬ 
tion |i). 

While our results indicate that these models all have interesting interpreta¬ 
tions as perturbed boundary conformal field theories, the checks performed so far 
should be considered as preliminary, and more work will be needed to put these 
results on a firmer footing. Two immediate open problems are fixing completely 
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the boundary UV/IR relations, and determining the boundary bound-state con¬ 
tent of each ADET model. The effort involved should be rewarded both by the 
well-established role in condensed matter physics taken by some of these theories, 
and by the mathematical insights into general integrable boundary theories their 
further study should bring. 
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A Diagonal gi x 01/02 coset data 

In general, to extract the g/f) coset conformal theory from the 0 WZW model, 
we need to decompose the representations, fi, of g into a direct sum of represen¬ 
tations, u, of f): 

/X 1-^ (A.l) 

U 

This decomposition corresponds to the character identity 

XmCt) = Xu{T)bfj,i,{T) (A.2) 

where x^^ and Xu are the characters of the g and f) representations and u 
respectively and the branching function, 6^jy(r), is the character of the coset 
theory. Let 11 denote the projection matrix giving the explicit projection of every 
weight of 0 onto a weight of f). Clearly for a coset character to be non zero we 
must have 

U/^o - izo G UQ (A.3) 

where Q is the root lattice of g. For our diagonal cosets, since ® Q) = Q 
this selection rule is particularly simple: 

y-o + ^0 — Po ^ Q (A-4) 
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where /Iq, i^o and pg are weights of 0. 

The group of outer automorphisms of 0, 0(0), permutes the fundamental 
weights in such a way as to leave the extended Dynkin diagram invariant. The 
action of an element, A G 0(0), on a modular S-matrix element of 0 at some level 
I is 

= (A.5) 

The modular S-matrix for a diagonal coset theory has the form 


c _ 0(1) 0(1) c( 2 )* 

Under the 0(0) action it transforms as 


^{Ap,Aiy,Ap}{fi' 


g27ri(AAo,/Xo-l-»^0-Po) g 




S 




(A.6) 


(A.7) 


since pg -|- p'q — pg G Q [57]. This suggests that 5 is a degenerate matrix, 
which cannot be the case since it represents a modular transformation. We must 
therefore remove this degeneracy by identifying the fields 


{Ap, Au; Ap} = {p, u; p} VA G 0(0). (A.8) 


For these diagonal cosets, the orbit of every element A of 0(0) has the same 
order, N, which is simply the order of the global symmetry group of the model. 
With this field identification, we find that every field in the theory appears with 
multiplicity N. To remedy this, the partition function must be divided by this 
multiplicity, which has the effect of introducing N into the coset modular S-matrix 
[58], 

(A,9) 

More information on field identifications can be found, for example, in [59]. 

For the level 1 representations of simply laced affine Lie algebras, the charac¬ 
ters have a particularly simple form [60, 61] : 

where q = exp(27rfr) and the generalised 9 and Dedekind rj functions are 

Opiq) = Y1 (A.ii) 

00 

v{q) = (A.12) 

n=l 




(A.icSj 
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where P and Q are the weight lattice and root lattice respectively [61, 62], and 
the r] function becomes 

vi--) = (A.14) 

r 

The modular S-matrix for level 1 is therefore 


S 


( 1 ) _ 

flU — 


1 

VWM 


g-2vri(i/o,Mo) 


(A.15) 


For the coset description of the (jf-function, eq.(7.20), 
which is simply 


c(l) _ 


1 

VWlQl' 


we need to compute S'q^ 
(A.16) 


It is important to note that P/Q is isomorphic to the centre of the group un¬ 
der consideration [63, 59]. The group of field identifications of the coset is also 
isomorphic to the centre of this group so consequently eq.(A.16) becomes 


o(l) _ J_ 

0 /^ ^ 


for all /X. 


(A.17) 


The (/-function can now be written in terms of the level 2 modular S matrices 
only; 




S, 


( 2 ) 

Op 


S. 


( 2 ) 

00 


(A.18) 


As we can see from (A.18), identifying the conformal g'-function value will 
only pin down the level 2 representation. For Eg, since there is only one possible 
level 1 representation (with Dynkin labels n* = 0, i = 1,..., 8), by specifying the 
level 2 representation we fix the coset field. However, for other cases, although the 
coset selection and identification rules (A.4),(A.8) do constrain the possible coset 
representations we are still left with some ambiguity in general. For example, 
for the Ar cases the selection and identification rules are quite simple and the 
possible cosets, given a fixed level 2 weight, are shown in table 6. The notation 
is as follows: Vj are level 1 weights and pj is the level 2 weight with Dynkin 

labels Uj = 1 for i = j and n* = 0, i = 1,..., r otherwise. The labels /Xg, uq and 
pQ now represent 0 with Dynkin labels n* = 0, i = 1,..., r. 

The coset fields for the A 2 (three-state Potts model) and (tricritical Ising 
model) cases, along with the corresponding boundary condition labels from [49] 
and [50] respectively, are given in tables 7 and 8 as concrete examples. 

It is useful to note that for the Ar, Dr, Eq and Ei models, 5oo = *S'op only 
when p = AO for some A G 0(0) [64] and for each such p there is a unique 
coset field, so the number of fields with (/-function equal to (/{y_t^jy^o} is equal to 
the size of the orbit of 0. On the other hand, E^ has no diagram symmetry, but 
it is also exceptional in that 5oo = Sop^ where P 2 has Dynkin labels n2 = 1, all 
other Hi = 0. Physically, this is to be expected as the two fields correspond to 
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Fixed level 2 weight 

Level 1 weight 

Level 1 weight 

Number of coset 

P 

P 

u 

fields {p, v; p} 

0 

0 

0 



Pi 

^ — 1, . . . r 

r + 1 

Pj 

Pi 

o' 

II 

T 


j = l,...[(r + l)/2] 


v>r+l-k, fc = 1,. ..r-j 

r + 1 


Table 6: Ar coset fields, indicating the number of distinct fields for each level 2 weight; 
in the first column, [x] denotes the integer part of x 


A 2 

Coset field 

Boundary label 

Level 2 weight label 


from [49] 

from table 2 

{[1,0,0], [1,0,0]; [2,0,0]} 

A 


{[0,0,1], [0,1,0]; [2, 0,0]} 

B 

0 

{[0,1,0], [0,0,1]; [2,0,0]} 

C 


{[1,0,0], [0,1,0]; [1,1,0]} 

AB 


{[0,0,1], [0,0,1]; [1,1,0]} 

BC 

Pi 

{[0,1,0], [1,0,0]; [1,1,0]} 

AC 



Table 7; A 2 coset fields with the corresponding boundary labels and level 2 weight labels 
from table 2. The weights are given in terms of Dynkin labels [no, ni, n2, •..] 


the two fixed boundary conditions (—) and (+), which clearly must have equal g- 
function values. (This exceptional equality is discussed from a more mathematical 
perspective in, for example, [64].) 

From (A.5) it is clear that 5'o(Ap) = Sop is also true for p / 0. For the A and 
E models at level 2, we find that these are the only cases where Sop. = Sop^, for 
the Dr models there is more degeneracy. 
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Ev 

Coset field 

Boundary label 

Level 2 weight label 


from [50] 

from table 2 

{[1,0 ,...,0],[1,0,...,0],[2,0,...,0]} 

(-) 

0 

{[0,1,0,...],[0,1,0,...];[2,0,...,0]} 

(+) 


{[1,0,...,0],[0,1,0,...];[!,1,0,...,0]} 

id) 

Pi 

{[1,0,...,0],[1,0,...,0];[0,0,1,0,...]} 

(-0) 

P2 

{[0,1,0,...], [0,1,0,...]; [0,0,1,0,...]} 

(0+) 


{[1,0,...,0], [0,1,0,...]; [0,0,0,1,0,...]} 

(0) 

P 3 


Table 8; Ey coset fields with the corresponding boundary labels and level 2 weight labels 
from table 2. The weights are given in terms of Dynkin labels [no, ni, n 2 ,...] 


References 

[1] Al.B. Zamolodchikov, talk at the 4th Bologna workshop on conformal and 
integrable models, July 3, 1999 [unpublished]. 

[2] P. Dorey, A. Pocklington, R. Tateo and G.M.T. Watts, ‘TEA and TCSA 
with boundaries and excited states’, Nucl. Phys. B 525 (1998) 641 
[arXiv:hep-th/9712197]. 

[3] Z. Bajnok, L. Palla and G. Takacs, ‘Spectrum and boundary energy in 
boundary sine-Gordon theory’, Nucl. Phys. B 622 (2002) 565 

[ar Xiv: hep-th/0108157]. 

[4] E. Gorrigan and A. Taormina, ‘Reflection factors and a two-parameter 
family of boundary bound states in the sinh-Gordon model’, J. Phys. A 33, 
8739 (2000) [arXiv;hep-th/0008237]. 

[5] M. Kormos and L. Palla, ‘Some semi-classical issues in boundary 
sine-Gordon model’, J. Phys. A 35, 5471 (2002) [arXiv:hep-th/0201230]. 

[ 6 ] Al.B. Zamolodchikov, ‘Thermodynamic Bethe ansatz in relativistic models. 
Scaling three state Potts and Lee-Yang models’, Nucl. Phys. B 342, 695 
(1990). 

[7] T.R. Klassen and E. Melzer, ‘Purely elastic scattering theories and their 
ultraviolet limits’, Nucl. Phys. B 338 (1990) 485. 

[ 8 ] I. Affleck and A.W.W. Ludwig, ‘Universal noninteger ‘ground state 
degeneracy’ in critical quantum systems’, Phys. Rev. Lett. 67 (1991) 161. 

[9] A. LeClair, G. Mussardo, H. Saleur and S. Skorik, ‘Boundary energy and 
boundary states in integrable quantum field theories’, Nucl. Phys. B 453 
(1995) 581 [arXiv;hep-th/9503227]. 


41 




[10] P. Dorey, I. Runkel, R. Tateo and G.M.T. Watts, ‘g-function flow in 
perturbed boundary conformal field theories’, Nucl. Phys. B 578 (2000) 85 
[arXiv:hep-th/9909216]. 

[11] P. Dorey, D. Fioravanti, C. Rim and R. Tateo, ‘Integrable quantum field 
theory with boundaries: the exact g-function’, Nucl. Phys. B 696 (2004) 445 
[arXiv:hep-th /0404014]. 

[12] H.W. Braden, E. Corrigan, P. Dorey and R. Sasaki, ‘Affine Toda field 
theory and exact S matrices’, Nucl. Phys. B 338 (1990) 689. 

[13] V.A. Fateev and A.B. Zamolodchikov, ‘Conformal field theory and purely 
elastic S matrices’, Int. J. Mod. Phys. A 5, 1025 (1990). 

[14] P. Dorey, ‘Root systems and purely elastic S matrices’, Nucl. Phys. B 358 
(1991) 654. 

[15] P. Dorey, ‘Root systems and purely elastic S matrices. 2’, Nucl. Phys. B 
374 (1992) 741 [arXiv:hep-th/9110058]. 

[16] A. Fring and D.I. Olive, ‘The fusing rule and the scattering matrix of affine 
Toda theory’, Nucl. Phys. B 379 (1992) 429. 

[17] F. Ravanini, R. Tateo and A. Valleriani, ‘Dynkin TBAs’, Int. J. Mod. Phys. 
A 8 (1993) 1707 [arXiv:hep-th/9207040]. 

[18] A.B. Zamolodchikov, ‘Integrals of motion and S matrix of the (scaled) T = 
T(C) Ising model with magnetic field’, Int. J. Mod. Phys. A 4, 4235 (1989). 

[19] A.B. Zamolodchikov, ‘Integrable field theory from conformal field theory’. 
Proceedings of Taniguchi Symposium, Kyoto (1988). 

[20] P. Dorey, ‘Exact S matrices ’, in the proceedings of the 1996 Eotvos 
Graduate School [arXiv:hep-th/9810026]. 

[21] F.A. Smirnov, ‘The perturbated c < 1 conformal field theories as reductions 
of sine-Gordon model’, Int. J. Mod. Phys. A 4 (1989) 4213. 

[22] J.L. Cardy and G. Mussardo, ‘S matrix of the Yang-Lee edge singularity in 
two-dimensions’, Phys. Lett. B 225 (1989) 275. 

[23] P.G.O. Freund, T.R. Klassen and E. Melzer, ‘S matrices for perturbations 
of certain conformal field theories’, Phys. Lett. B 229, 243 (1989). 

[24] P. Dorey and R. Tateo, ‘Excited states in some simple perturbed conformal 
field theories’, Nucl. Phys. B 515 (1998) 575 [arXiv;hep-th/9706140]. 

[25] A. Fring and R. Koberle, ‘Eactorized scattering in the presence of reflecting 
boundaries’, Nucl. Phys. B 421 (1994) 159 [arXiv:hep-th/9304141]. 

[26] S. Ghoshal and A.B. Zamolodchikov, ‘Boundary S matrix and boundary 
state in two-dimensional integrable quantum field theory’, Int. J. Mod. Phys. 
A 9 (1994) 3841 [Erratum-ibid. A 9 (1994) 4353] [arXiv:hep-th/9306002]. 


42 



[27] R. Sasaki, ‘Reflection bootstrap equations for Toda field theory’, in the 
proceedings of the conference, Interface between physics and mathematics, 
eds. W. Nahm and J.-M. Shen [arXiv;hep-th/9311027]. 

[28] P. Dorey, unpublished. 

[29] E. Corrigan, P. Dorey, R.H. Rietdijk and R. Sasaki, ‘Affine Toda field 
theory on a half line’, Phys. Lett. B 333 (1994) 83 [arXiv;hep-th/9404108]. 

[30] E. Corrigan, P. Dorey and R.H. Rietdijk, ‘Aspects of affine Toda field 
theory on a half line’. Prog. Theor. Phys. Suppl. 118 (1995) 143 
[arXiv:hep-th/9407148]. 

[31] V.A. Eateev, ‘Normalization factors, reflection amplitudes and integrable 
systems ’, [arXiv:hep-th/0103014]. 

[32] C. Zambon, York PhD Thesis 2004, unpublished. 

[33] P. Dorey and F. Ravanini, ‘Staircase models from affine Toda field theory’, 
Int. J. Mod. Phys. A 8 (1993) 873 [arXiv;hep-th/9206052]. 

[34] S. Ghoshal, ‘Bound state boundary S matrix of the sine-Gordon model’, Int. 
J. Mod. Phys. A 9 (1994) 4801 [arXiv;hep-th/9310188]. 

[35] H.J. de Vega and A. Gonzalez Ruiz, ‘Boundary K matrices for the six 
vertex and the n(2n-l) A(n-l) vertex models’, J. Phys. A 26 (1993) L519 
[arXiv:hep-th/9211114]. 

[36] J.L. Cardy, ‘Conformal invariance and surface critical behavior’, Nucl. 
Phys. B 240 (1984) 514. 

[37] Al.B. Zamolodchikov, ‘On the thermodynamic Bethe ansatz equations for 
reflectionless ADE scattering theories’, Phys. Lett. B 253 (1991) 391. 

[38] P. Fendley, H. Saleur and N.P. Warner, ‘Exact solution of a massless scalar 
field with a relevant boundary interaction’, Nucl. Phys. B 430, (1994) 577 
[ar Xiv: hep-th/9406125]. 

[39] V.V. Bazhanov, S.L. Lukyanov and A.B. Zamolodchikov, ‘Integrable 
structure of conformal field theory, quantum KdV theory and 
thermodynamic Bethe ansatz’, Commun. Math. Phys. 177, (1996) 381 
[ar Xiv: hep-th/9412229]. 

[40] V.V. Bazhanov, A.N. Hibberd and S.M. Khoroshkin, ‘Integrable structure 
of W(3) conformal field theory, quantum Boussinesq theory and boundary 
affine Toda theory’, Nucl. Phys. B 622, (2002)475 [arXiv;hep-th/0105177]. 

[41] F.A. Smirnov, ‘Reductions of the sine-Gordon model as a perturbation of 
minimal models of conformal field theory’, Nucl. Phys. B 337 (1990) 156. 

[42] R. Chatterjee, ‘Exact partition function and boundary state of 2-D massive 
Ising field theory with boundary magnetic field’, Nucl. Phys. B 468 (1996) 
439 [arXiv:hep-th/9509071]. 


43 



[43] V. Fateev, S.L. Lukyanov, A.B. Zamolodchikov and Al.B. Zamolodchikov, 
‘Expectation values of local fields in Bullough-Dodd model and integrable 
perturbed conformal field theories’, Nucl. Phys. B 516 (1998) 652 
[arXiv: hep-th/9709034]. 

[44] R. Guida and N. Magnoli, ‘Vacuum expectation values from a variational 
approach’, Phys. Lett. B 411 (1997) 127 [arXiv;hep-th/9706017]. 

[45] P. Di Francesco, P. Mathieu and D. Senechal, ‘Conformal field theory’. 
Springer, 1997. 

[46] P. Goddard, A. Kent and D.I. Olive, ‘Virasoro algebras and eoset spaee 
models’, Phys. Lett. B 152 (1985) 88. 

[47] T. Gannon, ‘Algorithms for affine Kae-Moody algebras’, 

[ar Xiv: hep-th/0106123]. 

[48] B. Schellekens, http://www.nikhef.nl/~t58/kac.html. 

[49] 1. Affleck, M. Oshikawa and H. Saleur, ‘Boundary critical phenomena in the 
three-state Potts model’, [arXiv:cond-mat/9804117]. 

[50] L. Ghim, ‘Boundary S-matrix for the trieritieal Ising model’, Int. J. Mod. 
Phys. A 11 (1996) 4491 [arXiv;hep-th/9510008]. 

[51] V.A. Fateev, ‘The exact relations between the coupling constants and the 
masses of particles for the integrable perturbed conformal field theories ’, 
Phys. Lett. B 324 (1994) 45. 

[52] S. Fredenhagen, ‘Organizing boundary RG flows’, Nucl. Phys. B 660 (2003) 
436 [arXiv:hep-th/0301229]. 

[53] 1. Affleck, ‘Edge critical behaviour of the 2-dimensional tri-critical Ising 
model’, J. Phys. A 33 (2000) 6473 [arXiv;cond-mat/0005286]. 

[54] K. Graham, ‘On perturbations of unitary minimal models by boundary 
condition changing operators’, JHEP 0203, 028 (2002) 

[arXiv:hep-th/0111205]. 

[55] J.L. Gardy, ‘Boundary conditions, fusion rules and the Verlinde formula’, 
Nucl. Phys. B 324 (1989) 581. 

[56] H. Saleur and M. Bauer, ‘On some relations between local height 
probabilities and conformal invariance’, Nucl. Phys. B 320 (1989) 591. 

[57] D. Gepner, ‘Field identification in coset conformal field theories’, Phys. 
Lett. B 222 (1989) 207. 

[58] A.N. Schellekens and S. Yankielowicz, ‘Field identification fixed points in 
the coset construction’, Nucl. Phys. B 334 (1990) 67. 

[59] G.R. Ahn and M.A. Walton, ‘Field identifications in coset conformal 
theories from projection matrices’, Phys. Rev. D 41 (1990) 2558. 


44 



[60] LB. Frenkel and V.G. Kac, ‘Basic representations of affine Lie algebras and 
dual resonanee models’, Invent. Math. 62 (1980) 23. 

[61] V.G. Kac and D.H. Peterson, ‘Infinite dimensional Lie algebras, theta 
functions and modular forms’, Adv. Math. 53 (1984) 125. 

[62] V.G. Kac, ‘Infinite dimensional Lie algebras’, Cambridge University Press, 
1990. 

[63] D. Bernard, ‘String charaeters from Kae-Moody automorphisms’, Nucl. 
Phys. B 288 (1987) 628. 

[64] J. Fuchs, ‘Simple WZW eurrents’, Commun. Math. Phys. 136 (1991) 345. 


45 



